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I.  STATEMENT  OF  THE  PROBLEM 


The  fundamental  problem  in  penetration  mechanics  may  be  stated  as  (Wright  and 
Frank1):  Given  a  projectile,  a  target,  and  details  of  their  initial  geometry,  kinematics, 
and  materials;  determine  whether  or  not  the  target  will  be  perforated  upon  impact.  If 
perforated,  find  the  residual  characteristics  of  the  projectile  and  the  target,  and  if  not, 
ascertain  the  shape  and  depth  of  the  hole.  Of  the  three  approaches,  namely,  simple  data 
correlation,  engineering  models,  and  numerical  simulation  used  to  analyze  a  penetration 
problem,  we  have  concentrated  on  the  third  approach  and  developed  robust  and  efficient 
algorithms  capable  of  analyzing  details  of  thermomechanical  deformations  of  the  penetrator 
and  the  target. 

Factors  that  play  a  significant  role  during  the  penetration  of  metal  targets  by  projectiles 
include  material  properties,  impact  velocity,  projectile  shape,  target  support  position,  and 
relative  dimensions  of  the  target  and  the  projectile.  Recently,  emphasis  has  been  placed  on 
kinetic  energy  penetrators,  which  for  terminal  ballistic  purposes  may  be  considered  as  long 
metal  rods  traveling  at  high  speeds.  For  impact  velocities  in  the  range  of  2  to  10  km/sec, 
compressible  hydrodynamic  flow  equations  can  be  used  to  describe  adequately  the  impact 
and  penetration  phenomena,  because  large  stresses  occurring  in  hypervelocity  impact  permit 
one  to  neglect  the  rigidity  and  compressibility  of  the  striking  bodies.  Models,  which  require 
the  use  of  the  Bernoulli  equation  or  its  modification  to  describe  this  hypervelocity  impact, 
have  been  proposed  by  Birkhoff  et  al.2  and  Pack  and  Evans.3  At  ordnance  velocities  (0.5  to 
2  km/sec),  material  strength  becomes  an  important  parameter.  Allen  and  Rogers4  modified 
the  Pack  and  Evans3  flow  model  by  representing  the  strength  as  a  resistive  pressure.  This  idea 
was  taken  further  by  Alekseevskii5  and  Tate,6,7  who  considered  separate  resistive  pressures 
for  the  penetrator  and  the  target.  These  resistive  pressures  are  empirically  determined 
quantities,  and  the  predicted  results  depend  strongly  upon  the  assumed  values  of  these 
pressures.  As  described  lucidly  by  Wright8  in  his  survey  article  on  long  rod  penetrators, 
Tate’s  model  is  difficult  to  use  for  quantitative  purposes,  because  the  strength  parameters 
depend  upon  the  velocity  of  impact  and  the  particular  combination  of  materials  involved. 

The  one-dimensional  theories  ignore  the  lateral  motion,  plastic  flow,  and  the  detailed 
dynamic  effects.  The  paper  by  Backman  and  Goldsmith9  is  an  authoritative  review  of 
the  open  literature  on  ballistic  penetration,  containing  27S  reference  citations  from  the 
1800’s  to  1977.  They  describe  different  physical  mechanisms  involved  in  the  penetration 
and  perforation  processes,  and  also  discuss  a  number  of  engineering  models.  Jonas  and 
Zukas10  reviewed  various  analytical  methods  for  the  study  of  kinetic  energy  projectile- 
armor  interaction  at  ordnance  velocities  and  placed  particular  emphasis  on  three-dimensional 
numerical  simulation  of  perforation.  Anderson  and  Bodner11  have  recently  reviewed 
engineering  models  for  penetration  and  some  of  the  major  advances  in  hydrocode  modeling 
of  penetration  problems.  Three  books,12-14  published  during  the  past  few  years,  include 
extensive  discussions  of  the  engineering  models,  experimental  techniques,  and  analytical 
modeling  of  ballistic  perforation. 

Manganello  and  Abbot,15  Wingrove  and  Wulf,16  and  Recht17  observed  that  the 
penetration  resistance  of  some  armor  materials  is  reduced,  even  though  these  materials 
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mechanics  formulated  in  either  Lagrangian  (HEMP)  or  Eulerian  (HELP)  coordinates.  In 
the  HEMP  code,  the  finite  difference  operators  are  centered  in  space  and  time  to  give  second 
order  accuracy.  The  von  Mises  yield  criterion  and  the  incremental  theory  of  plasticity  are 
used  to  derive  the  incremental  stress-strain  relations.  The  material  model  employed  in  the 
HELP  code  includes  the  Tillotson  equation  of  state  modified  to  give  a  smooth  transition 
between  condensed  and  expanded  states,  the  von  Mises  yield  criterion  with  yield  strength 
taken  to  increase  with  hydrostatic  pressure  and  decrease  linearly  with  the  increase  in  the 
temperature,  and  a  failure  criterion  based  upon  either  a  minimum  value  of  the  mass  density 
or  a  critical  value  of  the  principal  stress  in  tension. 

Batra  and  Chen35  recently  developed  a  simple  technique  with  which  to  analyze 
approximately  the  axisymmetric  steady  state  deformations  of  a  rigid/viscoplastic  target 
being  penetrated  by  a  long  rigid  cylindrical  rod  with  a  hemispherical  nose.  They  assume  a 
kinematically  admissible  velocity  field  that  satisfies  the  balance  of  mass,  all  of  the  essential 
boundary  conditions,  and  traction  boundary  conditions  on  the  axis  of  symmetry  and  on  the 
target/penetrator  interface,  and  they  evaluate  various  parameters  in  the  presumed  velocity 
field  by  minimizing  the  error  in  the  satisfaction  of  the  balance  of  linear  momentum.  The 
computed  results  reveal  that  the  leading  term  in  the  proposed  velocity  field  gives  a  good 
solution  that  is  reasonably  close  to  the  previously  computed  finite  element  solution  of  the 
problem.  Batra  and  Gobinath36  analyzed  the  steady  state  axisymmetric  penetration  problem 
for  a  deformable  penetrator  and  target,  and  modeled  the  material  for  each  as  thermally 
softening,  but  strain  and  strain-rate  hardening.  They  found  that  the  bottom  part  of  the 
target/penetrator  interface  is  ellipsoidal  rather  than  hemispherical.  The  peak  pressures  in 
the  penetrator  near  the  stagnation  point  approach  4.58crop  and  that  in  the  target  I4cr0t  when 
Vo-p/vot  —  3-06.  Here  crop  and  a0t  equal,  respectively,  the  yield  stress  in  a  quasistatic  simple 
compression  test  for  the  penetrator  and  target  materials.  The  axial  resisting  force  on  the 
penetrator  equaled  8.91  F,  11.52  F,  and  14.51  F  (F  =  nr20aop )  for  stagnation  point  speeds 
of  450  m/s,  500  m/s,  and  550  m/s,  respectively.  A  significant  contribution  to  the  resisting 
force  is  made  by  the  consideration  of  strain-rate  hardening  effects.  When  the  penetrator  and 
target  materials  are  modeled  as  rigid/perfectly  plastic,  Gobinath  and  Batra37  found  that  the 
resistive  pressure  terms  in  the  modified  Bernoulli’s  equation  depend  upon  the  ratio  of  the 
mass  densities  of  the  target  and  the  penetrator  as  well  as  on  the  penetration  speed. 

The  analysis  of  the  steady  state  axisymmetric  deformations  of  an  elastic/perfectly 
plastic  target  being  penetrated  by  a  fast  moving  rigid  cylindrical  rod  by  Jayachandran  and 
Batra38  has  revealed  that  the  consideration  of  elastic  effects  reduces  the  value  of  the  peak 
hydrostatic  pressure  acting  at  the  stagnation  point,  the  axial  resisting  force  experienced  by 
the  penetrator,  and  the  target  resistance  parameter  appearing  in  the  modified  Bernoulli 
equation  nearly  28%,  25%,  and  25%,  respectively.  Batra  and  Jayachandran39  modeled 
the  thermoviscoplastic  response  of  the  target  material  by  flow  rules  due  to  Litonski-Batra, 
Bodner-Partom,40  and  Brown,  Kim,  and  Anand.41  They  calibrated  the  three  flow  rules 
against  a  hypothetical  plane  strain  compression  test  performed  at  a  nominal  strain-rate  of 
3,300  sec-1  and  then  studied  the  axisymmetric  penetration  problem.  They  found  that  the 
three  constitutive  relations  gave  nearly  the  same  value  of  the  resisting  force  acting  on  the 
penetrator,  temperature  rise  of  material  particles  in  the  vicinity  of  the  target/penetrator 
interface,  and  other  macroscopic  measures  of  deformation,  such  as  the  effective  stress  and 
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Fig.  2).  This  layer  adjoins  the  surface  of  the  hole  punched  in  the  plate  and  piastic 
deformations  of  the  material  within  the  layer  are  very  intense.  Even  though  the  mesh  has 
been  refined  adaptively,  the  scarcity  of  computational  resources  has  limited  our  ability  to 

decipher  the  width  of  the  shear  band. 

In  an  attempt  to  understand  the  susceptibility  of  tungsten  and  depleted  uranium 
penetrators  to  adiabatic  shear  banding,  we  have  investigated  the  development  of  shear  bands 
in  a  rectangular  block  of  these  materials  deformed  in  plane  strain  compression  at  a  nominal 
strain-rate  of  5000/s.  It  is  assumed  that  the  deformations  are  symmetric  about  the  horizontal 
and  vertical  centroidal  axes;  therefore,  deformations  of  a  quarter  of  the  domain  are  studied. 
The  undeformed  region  is  divided  into  104  uniform  rectangular  elements  and  it  is  assumed 
that  the  yield  strength  of  300  randomly  located  elements  in  the  first  quadrant  is  5%  lower 
than  that  of  the  remaining  elements.  In  each  case  the  material  response  is  modeled  by 
the  Johnson-Cook  relation  with  the  values  of  material  parameters  taken  from  Rajendran’s 
report.46  As  depicted  in  Figure  3,  the  pattern  of  shear  bands  formed  in  the  two  materials 
is  quite  different;  similar  patterns  were  obtained  for  100  and  200  randomly  located  weak 
elements.  In  each  case,47  shear  bands  initiated  in  tungsten  at  a  lower  value  of  the  average 
strain  than  that  in  depleted  uranium.  Batra,  Zhang,  and  Wright48  obtained  a  similar  result 
for  these  two  materials  deformed  in  torsion  which  agrees  with  the  experimental  observations 
of  Johnson  et  a/..49 

We  have  incorporated  the  element  erosion  algorithm  in  conjunction  with  the  Johnson- 
Cook  model  in  DYNA2D.  With  this  there  is  no  need  to  refine  the  mesh  adaptively  since 
severely  deformed  elements  are  eliminated  from  the  computations.  Computations  with 
identical  depleted  uranium  and  tungsten  penetrators  striking  at  normal  incidence  rolled 
homogeneous  armor  (RHA)  steel  plates  reveal  that  the  residual  speed  of  depleted  uranium 
penetrators  is  more  than  that  of  tungsten  penetrators  (cf.  Fig.  4).  This  suggests  that  the 
depleted  uranium  penetrators  should  penetrate  deeper  into  the  thick  target  than  identical 
tungsten  penetrators. 

We°°  have  also  proposed  a  four-stage  engineering  model  of  the  penetration/perforation 
by  hemispherical-nosed  rigid  cylindrical  rods  into  targets  whose  material  exhibits  strain- 
hardening  effects.  During  each  stage  of  the  penetration  process,  a  kinematically  admissible 
velocity  field  involving  one  or  more  unknown  parameters  is  assumed.  These  parameters 
are  determined  by  minimizing  the  rate  of  plastic  dissipation  from  the  velocity  field,  the 
incremental  deformations  of  the  target,  the  penetration  depth,  the  resisting  force  acting  on 
the  penetrator  and  hence  its  deceleration  are  evaluated.  A  simple  criterion  for  the  formation 
and  ejection  of  the  cylindrical  plug  in  the  target  has  been  proposed.  Computed  results  (cf. 
Fig.  5)  for  the  exit  speed  of  the  penetrator  and  when  the  targets  are  not  perforated  of  the 
penetration  depth  are  found  to  match  well  with  the  corresponding  test  values. 
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Fig.  2  Distribution  of  axial  velocity  on  a  horizontal  line  in  the  plate  during  its  perforation  by 
a  rigid  punch.  The  length  is  nondimensionalized  with  respect  to  the  penetrator  radius. 
H  =  distance  from  the  punch  face;  the  nondimensional  time  equals  the  nondimensional 
depth  of  the  penetrator  into  the  plate. 
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Fig.  4  Penetration  depth  vs.  time  for  depleted  uranium  and  tungsten  rods  striking  at  normal 
incidence  RHA  steel  plates. 
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Abstract — We  study  axisvmmetnc  thermomechanical  deformations  of  a  thick  target  being  penetrated 
by  a  fast-moving  rigid  cylindrical  rod  with  a  hemispherical  nose,  and  presume  that  target  deformations 
appear  steady  to  an  observer  situated  at  the  penetrator  nose  tip.  Both  isotropic  and  kinematic 
hardening  of  the  target  material  are  considered.  It  is  found  that  kinematic  hardening  increases  the 
normal  stress  acting  on  the  penetrator  nose  surface  and  the  temperature  rise  of  target  particles 
abutting  the  penetrator.  However,  the  value  of  the  hydrostatic  pressure  at  a  point  in  the  deforming 
target  region  is  affected  very  little  by  the  consideration  of  kinematic  hardening.  For  suitable  values  of 
material  parameters  appearing  in  the  evolution  equation  of  the  back-stress,  the  computed  values  of 
the  back-stress  at  target  particles  abutting  the  penetrator  nose  surface  equal  three  times  the  yield 
stress  of  the  target  material  in  a  quasistatic  simple  compression  test. 


1.  INTRODUCTION 

During  the  penetration  of  a  thick  target  by  a  fast-moving  cylindrical  rod,  target  and  penetrator 
material  particles  in  the  vicinity  of  the  target/penetrator  interface  are  deformed  severely,  and 
are  also  heated  up  significantly.  Consequently,  the  material  undergoes  microstructural  changes 
such  as  the  generation/annihilation  of  dislocations,  dynamic  recovery  and  recrystallization, 
development  of  texture,  nucleation  and  growth  of  microcracks  and  voids,  and  possibly  the 
development  of  shear  bands  that  form  during  the  intense  plastic  deformations  of  a  material, 
especially  at  high  strain  rates.  One  way  to  account  for  these  microstructural  changes  is  to  use 
constitutive  equations  which  employ  a  suitable  number  of  scalar  and  tensor  valued  internal 
variables  (e.g.  see  Coleman  and  Gurtin  [1],  Chan  et  al.  [2],  Inoue  [3],  Lubliner  [4],  and  Anand 
[5]).  Here  we  use  one  scalar  variable  to  describe  the  isotropic  hardening,  and  a  traceless 
symmetric  second-order  tensor,  also  known  as  the  back  stress  tensor,  to  account  for  the 
kinematic  hardening  of  the  material.  The  Litonski-Batra  constitutive  relation  (e.g.  see  Batra 
and  Jayachandran  [6])  is  modified  to  incorporate  the  kinematic  hardening  and  used  herein  to 
study  the  thermomechanical  axisymmetric  deformations  of  the  target.  The  hemispherical  nosed 
cylindrical  penetrator  is  assumed  to  be  rigid  and  the  target  deformations  steady,  as  seen  by  an 
observer  situated  at  the  penetrator  nose  tip  and  moving  with  it.  We  note  that  Batra  and 
Jayachandran  [6]  recently  analyzed  thermomechanical  deformations  of  a  target  by  using  three 
constitutive  relations,  namely,  those  due  to  Litonski-Batra,  Bodner-Partom  [7],  and  Brown  et 
al.  [8],  Each  of  these  was  calibrated  to  give  almost  identical  effective  stress  vs  logarithmic  strain 
curves  for  a  block  made  of  target  material  and  deformed  in  plane  strain  compression  at  an 
average  strain-rate  of  3300  s-1.  Even  though  these  constitutive  relations  account  for  the 
evolution  of  the  microstructural  changes  in  different  ways,  they  gave  essentially  identical  results 
for  the  resisting  force  experienced  by  the  penetrator,  normal  stress  on  the  penetrator  nose 
surface,  and  the  distribution  of  the  tangential  speed  and  the  second-invariant  of  the  strain-rate 
tensor  on  the  penetrator  nose  surface. 

This  work  is  in  the  spirit  of  the  one  initiated  by  Batra  and  Wright  [9],  and  is  aimed  at 
providing  guidelines  for  selecting  and  improving  upon  the  previously  used  kinematically 
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Here  v  is  the  velocity  of  a  material  particle,  a  the  Cauchy  stress  tensor  at  the  present  location 
of  a  material  particle,  p  the  mass  density,  q  the  heat  flux,  D  the  stretching  tensor,  and  W  the 
spin  tensor.  The  balance  laws  (1),  (2),  and  (3)  are  written  in  the  Eulerian  description  of 
motion,  and  the  operators  grad  and  div  denote  the  gradient  and  divergence  operators  on  fields 
defined  in  the  present  configuration.  Equation  (1)  implies  that  both  the  elastic  and  plastic 
deformations  of  the  target  are  assumed  to  be  isochoric.  Equation  (5)  is  the  Fourier  law  of  heat 
conduction  with  k  the  thermal  conductivity  and  9  the  temperature  rise  of  a  material  particle. 
Equation  (6)  is  the  presumed  constitutive  relation  for  the  specific  internal  energy  U,  wherein  c 
is  the  specific  heat.  In  equations  (7)— ( 13),  a  is  the  Cauchy  stress  tensor,  s  its  deviatoric  part,  p 
the  hydrostatic  pressure  not  determined  by  the  deformation  history,  B  the  traceless  symmetric 
second-order  tensor  used  to  account  for  the  kinematic  hardening  of  the  material,  and  ip  is  a 
scalar  internal  variable  that  accounts  for  the  isotropic  hardening  of  the  material.  The  evolution 
equation  (10)  for  B  has  been  proposed  by  White  et  al.  [43],  Equation  (9)  with  p  given  by 
equation  (11)  is  a  generalization  of  the  Litonski-Batra  law  to  account  for  the  kinematic 
hardening  of  the  material.  Here  a„  is  the  yield  stress  in  a  quasistatic  simple  tension  or 
compression  test,  parameters  b  and  m  characterize  the  strain-rate  sensitivity  of  the  material,  v 
its  thermal  softening,  and  b  its  workhardening.  In  a  quasistatic  simple  tension  or  compression 
test, 


\  xpj 

describes  the  stress-strain  curve,  where  ip  is  now  interpreted  as  the  plastic  strain.  In  a  dynamic 
test,  the  effect  of  the  history  of  deformation  upon  the  present  state  of  deformation  is  accounted 
for  through  the  parameter  ip. 

We  nondimensionalize  variables  by  scaling  stress-like  quantities  by  cr„,  length  by  rn,  time  by 
(r0/u„),  and  the  temperature  by  the  reference  temperature  9r,  defined  by 

d'^ojpc.  (14) 


Here  r„  equals  the  radius  of  the  cylindrical  part  of  the  penetrator  and  v„  the  penetration  speed. 
In  terms  of  nondimensional  variables,  the  aforestated  governing  equations  become 


where 


divv  =  0,  (15.1) 

-  gradp  +  div  s  =  a{v  •  grad)v.  (15.2) 

s  -  B  +  Py((?  ■  grad)s  +  sW  -  Ws)  =  2(3D,  (15.3) 

tr(oDp)  +  <5  div(grad  9)  =  (v  •  grad)0,  (15.4) 

(v  ■  grad)B  +  BW-WB  +  ?2/B  =  s,Dp,  (15.5) 

(v  •  grad) xp  =  tr(oDp)/(  1  +— )  ,  (15.6) 

/  \  ip{)/ 


pry,  uv„ 

a  - - ,  p  = - . 

On  cr,,r0 


Y  —  —  ■  and  6  = - - 

G  pcv„r„ 


(15.7) 


are  nondimensional  numbers.  Henceforth,  we  will  use  nondimensional  variables  only.  For  a 
given  problem,  a,  y,  and  d  are  constants,  but  p  varies  from  point  to  point  in  the  deforming 
region  because  of  the  variation  in  p.  The  value  of  a  signifies  the  importance  of  inertia  forces 
relative  to  the  flow  stress  of  the  material,  and  may  be  thought  of  as  the  reciprocal  of  the 
Reynolds  number  in  a  viscous  fluid.  The  values  of  y  and  6  give  the  effect  of  material  elasticity 
and  heat  conduction,  respectively.  For  typical  penetration  problems  involving  long  rod 
penetrators,  <5  is  of  the  order  of  10~5;  hence  target  deformations  may  be  considered  adiabatic. 
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—  =  0  on  the  surface  AB,  (16.4) 

V>=0,  p  =0,  sn  =  0,  J9e  =  0,  szz  =  0, 

Bzz  =  0,  Btz  =  0  on  the  bounding  surface  EFA, 

(16.5) 
39 

2r  ~  0  on  the  axis  of  symmetry  DE.  (16.6) 

Here  n  and  t  denote,  respectively,  a  unit  normal  and  a  unit  tangent  vector  to  the  surface,  93  is 
an  average  temperature  of  the  penetrator,  hc  is  the  heat  transfer  coefficient  between  the 
penetrator  and  the  target,  and  T,  denotes  the  target/penetrator  interface.  The  boundary 
conditions  (16)  incorporate  the  assumptions  that  T,  is  smooth,  there  is  no  interpenetration  of 
the  target  material  into  the  penetrator  and  vice  versa,  the  deformations  are  axisymmetric,  and 
the  bounding  surfaces  AB  and  EFA  are  far  removed  from  the  penetrator  nose  surface.  The 
boundary  conditions  (16.5)  on  p  and  components  of  s  and  B  on  the  surface  EFA  are  needed, 
since  we  solve  equations  (15.3)  and  (15.5)  for  s  and  B,  along  with  the  other  equations  for  p,  v, 
9,  and  rp;  e.g.  see  Shimazaki  and  Thompson  [44j. 

We  refer  the  reader  to  [10]  for  details  of  obtaining  a  finite  element  solution  of  the  problem, 
and  ensuring  that  the  region  R  studied  herein  is  adequate.  The  computer  code  used  to  analyze 
the  thermomechanical  problem  discussed  in  [6]  was  modified  to  include  the  effect  of  kinematic 
hardening. 


vr  =  0,  vz  —  —  1 ,  9  =  0, 

srz  =0,  B  „  =  0,  B  gg  =  0, 

arz  =0,  ur  =  0, 


e  highly 
ek  their 
■  a  finite 
m  it  the 

(16.1) 

(16.2) 

(16.3) 


3.  NUMERICAL  RESULTS 

In  an  attempt  to  study  the  effect  of  kinematic  hardening  on  the  solution  variables,  the  values 
of  material  parameters  $|  and  g2  *n  equation  (15.5)  were  varied  over  a  wide  range.  However, 
the  other  material  and  geometric  parameters  were  assigned  the  following  values,  taken  from 
[6],  for  an  HY-100  steel. 

p  =  7860  kg/m3,  ct„  =  405  MPa,  G  =  80GPa,  c  =  473J/kg°C, 

k  =  50  W/nTC,  /t  =  20W/m2oC,  0a  =  0,  ro=10mm,  6  =  10  s, 

v  =  1.2  x  10_3/°C,  xp()  =  0.1,  m  =  0.01,  n=0. 13,  a  =  2.0.  (17) 

Thus,  the  reference  temperature  used  to  nondimensionalize  the  temperature  rise  equals 
108. 9°C. 

Wang  and  Batra  [45]  have  recently  studied  the  initiation  and  growth  of  shear  bands  in  a 
thermally  softening  viscoplasitc  block  obeying  constitutive  relations  similar  to  equations 
(8)— (12)  and  deformed  in  plane  strain  compression  at  an  average  strain-rate  of  5000  s-1.  Their 
load-displacement  curves  (cf.  Fig.  10  of  [45])  for  the  homogeneous  block  show  that  an  increase 
in  the  value  of  §,  hardens  the  material  and  an  increase  in  the  value  of  £2  softens  it  in  the  sense 
that  the  load  required  to  compress  the  block  by  a  certain  amount  is  more  for  higher  values  of 

and  less  for  larger  values  of  |2.  Earlier  computations  by  Batra  and  Jayachandran  [6]  and 
Jayachandran  and  Batra  [19]  suggest  that  a  change  in  the  values  of  and  |2  should  affect  the 
deformations  of  the  target  in  an  analogous  manner. 

All  of  the  results  presented  below  and  values  of  variables  indicated  in  figures,  unless  stated 
otherwise,  are  nondimensional. 

3. 1  Results  with  %2  varied 

Figure  2  depicts  the  distribution  of  the  normal  stress,  temperature  rise,  tangential  speed,  and 
the  second  invariant  of  the  strain-rate  tensor,  also  referred  to  as  the  strain-rate  measure,  on  the 
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Fig.  3.  Distribution  of  (tr(B  Br))''2  upon  the  penetrator  nose  surface  for  =  1.0  and  four  different 
values  of  |2.  See  Fig.  2  for  the  legend  to  curves.  The  curve  for  §2  =  100  essentially  coincides  with  the 

horizontal  axis. 


for  the  hemispherical  nose  surface  was  found  to  equal  7.87  for  the  case  of  no  kinematic 
hardening,  and  9.61,  8.57,  7.86,  and  7.77  for  £2  =  0,  1,  10,  and  100,  respectively.  The 
dimensional  values  of  F  are  obtained  by  multiplying  the  nondimensional  ones  by  Note 

that  in  our  work  the  speed  of  penetration  is  kept  fixed,  and  whatever  additional  energy  is 
required  for  the  penetration  process  is  presumed  to  be  available.  The  noticeable  increase  in  the 
temperature  rise  for  §2  =  0,  1,  and  10  is  due  to  the  significant  values  of  B,  see  Fig.  3,  and  the 
observation,  verified  by  the  computed  results  that  tr(BDp)2;0.  Since  values  of  (s-B)  affect 
the  plastic  deformations  of  the  material,  an  increase  in  the  value  of  B  will  necessitate  higher 
values  of  s,  which  will  cause  more  plastic  working  and,  hence,  greater  temperature  rise.  The 
variation  of  the  strain-rate  measure,  axial  stress,  temperature  rise,  and  the  axial  velocity  on  the 
central  line  plotted  in  Fig.  4  reveals  that  the  consideration  of  kinematic  hardening  affects  these 
variables  at  points  situated  at  most  one  penetrator  radius  from  the  penetrator  nose  tip. 

An  integration  of  equation  (15.2)  along  the  central  streamline  (r  =  0)  gives 

^  av2  +  p  -  s..  -  2  f  — ^  dz  =  -a..(0).  (19) 

2  ~  Jo  dr 

Setting  z  =  0  and  comparing  the  result  with  Tate’s  equation  (33,  34],  we  get 

=  (20) 

where  Rt  equals  the  strength  parameter  for  the  target  in  Tate's  equation,  and  <f.z  is  the  value  of 
azz  at  the  stagnation  point.  The  computed  values  of  Rt  for  if,  =  0,  1,  10,  and  100  were  found  to 
be  10.63,  8.91,  8.05,  and  7.97,  respectively.  For  the  case  of  no  kinematic  hardening,  Rt  =  8.01. 
According  to  Tate  [33,  34], 


where  £  equals  Young’s  modulus  for  the  target  material.  Thus,  Tate’s  formula  gives  £,  =  6.64. 
We  note  that  for  f2  =  0  and  1,  the  values  of  sszz  were  higher  than  those  for  £2  —  10  and  100. 
Since  the  plastic  deformation  is  governed  by  (szz  -  Bzz),  higher  values  of  Bzz  necessitate  a 
corresponding  increase  in  the  values  of  szz.  For  §2  =  0,  1,  10,  and  100,  the  values  of  -jzz(0) 
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Fig.  6.  Distribution  of  (tr(B  Br))1 2  upon  the  penetrator  nose  surface  for  =  104  and  four  different 
values  of  |2.  See  Fig.  5  for  the  legend  to  curves.  The  curve  for  §2  =  108  coincides  with  the  horizontal 

axis. 
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Fig.  8.  Distribution  of  (tr(B  BT))1'2  upon  the  penetrator  nose  surface  for  §2  =  100  and  four  different 
values  of  See  Fig.  7  for  the  legend  to  curves.  The  curve  for  =0. 1  coincides  with  the  horizontal 

axis. 
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Fig.  9.  Distribution  of  the  normal  stress,  tangential  speed,  temperature  rise,  and  the  second  invariant 
of  the  strain-rate  tensor  on  the  hemispherical  nose  surface  of  the  penetrator  for  5,  =  106  and  four 

different  values  of  - no  kinematic  hardening; - §,  =  102;  . 2-,  =  10;  -•-§,  =  10s; 

- =  1.5  x10s. 
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Fig.  11.  Distribution  of  the  normal  stress,  tangential  speed,  temperature  rise,  and  the  second 
invariant  of  the  strain-rate  tensor  on  the  hemispherical  nose  surface  of  the  penetrator  for  ?,  =  1, 
,1:  =  100,  and  three  different  values  of  or; - a  =  2; - a  =  6;  and - a  =  10. 


from  the  crater  surface,  and  then  quite  slowly.  The  contours  of  IB  for  other  values  of  f ,  and  £2 
are  similar  to  those  shown  in  Fig.  10  and  are  not  included  herein. 

3. 3  Results  with  a  varied 

We  have  plotted  in  Fig.  11  the  distribution  on  the  penetrator  nose  surface  of  the  normal 
stress,  temperature  rise,  tangential  speed,  and  the  strain-rate  measure  for  a  =  2,  6,  and  10  with 

=  1  and  §2  =  100  kept  fixed.  The  normal  stress  distribution  on  the  penetrator  nose  surface 
resembles  that  computed  by  Batra  and  Wright  [9]  for  the  case  of  no  kinematic  hardening  in  the 
sense  that  its  value  at  the  point  for  which  \p  ~  46°  is  unaffected  by  the  value  of  a,  and  it 
increases  with  a  for  xp  <  46°  and  decreases  with  a  for  i p  >  46°.  We  note  that  Batra  and  Wright 
used  a  coarse  finite  element  mesh  consisting  of  6-noded  triangular  elements  and  considered  a 
smaller  target  region  than  that  studied  herein.  The  values  of  the  strain-rate  measure  increase 
with  or  for  ip  <  60°,  those  of  the  tangential  speed  increase  with  a  for  10°  <  a  <  80°,  and  the 
temperature  distribution  on  the  penetrator  nose  surface  is  affected  very  little  when  a  is 
increased  from  two  to  ten.  The  distribution  of  (tr(B  BT))1'2  on  the  penetrator  nose  surface, 
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stayed  undeformed  for  all  practical  purposes.  Here  we  take  the  penetrator  to  be  rigid  and  give 
details  of  the  deformation  fields  such  as  history  of  the  temperature  at  a  material  point,  history  of 
the  penetrator  speed,  and  the  distribution  of  the  temperature  and  velocity  field  in  the  deforming 
target  region.  For  the  sixteen  tests  simulated  herein,  the  computed  depth  of  penetration  has  been 
found  to  match  very  well  with  that  reported  by  Forrestal  et  al.  [5], 

We  refer  the  reader  to  the  review  articles  by  Backman  and  Goldsmith  [6],  Wright  and  Frank 
[7],  Anderson  and  Bodner  [8],  and  the  books  by  Zukas  et  al.  [9],  Blazynski  [10],  MaCauley  [11], 
and  Zukas  [12]  for  a  summary  of  the  available  literature  on  ballistic  penetration.  Awerbuch  and 
Bodner  [13],  Ravid  and  Bodner  [14],  Ravid  et  al.  [15],  Forrestal  et  al.  [16],  and  Batra  and  Chen  [17] 
have  proposed  engineering  models  of  different  complexity  in  which  a  rigid  penetrator  impacts 
a  deformable  terget.  Batra  [18]  has  pointed  out  that  the  target  deformations  during  the  steady 
state  phase  of  the  penetration  process  are  influenced  strongly  by  the  penetrator  nose  shape.  He 
assumed  the  target/penetrator  interface  to  be  smooth.  Here  we  also  investigate  the  effect  of  the 
nose  shape  on  transient  deformations  of  the  target  and  account  for  the  effect  of  the  frictional  force 
at  the  interface.  The  shapes  of  the  intensely  deforming  target  material  ahead  of  the  penetrator 
nose  for  ellipsoidal  and  blunt  nosed  penetrators  are  quite  different.  In  each  case  the  plastic 
deformations  of  the  target  material  are  large  enough  to  melt  a  thin  layer  of  the  target  material 
adjoining  the  target/penetrator  interface. 


2  Formulation  of  the  problem 

In  order  to  keep  the  paper  self-contained  we  summarize  below  the  pertinent  equations.  We 
presume  that  target  deformations  are  axisymmetric.  and  use  a  fixed  set  of  cylindrical  coordinates 
w'ith  r-axis  coincident  with  the  axis  of  symmetry  and  pointing  into  the  target  and  the  origin  at  the 
top  surface  of  the  undeformed  target  to  describe  its  deformations.  The  balance  laws  written  in  the 
referential  description  of  motion  and  governing  the  thermomechanical  deformations  of  the 
target  are 


(oi)  =  0,  (1.1) 

<j0v  =  Div  T.  (1,2) 

i?o<?  =  -Div  Q  +  tr  (TFT),  (1.3) 

where 

J  =  det  F,  F=  Gradx.  (2) 


x  is  the  present  position  of  a  material  particle  that  occupied  place  A'  in  the  reference 
configuration,  g  its  present  mass  density.  g0  its  mass  density  in  the  reference  configuration,  v  the 
present  velocity  of  a  material  particle.  T  the  first  Piola-Kirchhoff  stress  tensor,  e  the  specific 
internal  energy.  Q  the  heat  flux  per  unit  reference  area,  a  superimposed  dot  indicates  the  material 
time  derivative,  and  operators  Grad  and  Div  signify  the  gradient  and  divergence  of  field 
quantities  defined  in  the  reference  configuration.  The  balance  law's  1 1. 1)  —  ( 1.3)  are  supplemented 
by  the  following  constitutive  relations: 

<r=-p(o)l  +2 jiD,  T  =  —  o[F~  !)r,  (3.1,2) 
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Fig.  1.  A  schematic  sketch  of  the 
problem  studied 


.--direction,  and  we  reckon  time  from  the  instant  it  impacts  the  top  surface  of  the  target  at  normal 
incidence. 

In  order  to  solve  the  problem  numerically  we  analyze  deformations  of  the  finite  target  region 
ABCDEFA  shown  in  Fig.  1  and  impose  the  following  boundary  conditions  on  the  bounding 
surfaces: 

an  =  0.  q  ■  n  =  h(9  —  0J  on  FED.  (5.1.2) 

ar:  =  0,  vr  =  0.  qr  =  0  on  the  axis  of  symmetry  AB.  (5.3  -  5l 

q  n  =  0,  an  =  0  at  a  point  on  AF  where  the  penetrator  surface  is  not  in  contact 

with  the  deforming  target  region.  (5.6.  7) 

q  n  =  0,  [v ■  /i]  =  0 ,  e,  an  =  at  points  on  AF  where  the  penetrator  and  target 

surfaces  are  in  contact  with  each  other.  (5.8—  10) 

f„  =  —gVprn,  f  =  —  oF'e,.  q  n  =  0.  on  the  bounding  surface  BCD.  (5.1 1  -  13) 

Here  Vp  and  Vs  are  the  speeds  in  the  target  material  of  the  p- waves  and  s-waves.  respectively. 
i'„  and  v,  are  the  normal  and  tangential  components  of  the  velocity  of  a  target  particle  on  the 
bounding  surface  BCD.  and  /„  and  f,  equal,  respectively,  the  normal  and  tangential  forces  exerted 
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values  to  other  material  parameters: 

go  =  2710kg/m3,  ,4  =  68.9  GPa.  B  =  0.  k  =  120  Wm  "  1  C"  1 . 


:  875 Jkg" 1  C~l, 


22  C.  h  —  20 Wm"’  C 


p  =  1.5. 


According  to  our  constitutive  relations  (3.2)  and  (3. 3),  when  6  =  l  i;  =  melting  temperature  of  the 
material,  u  =  0.  the  material  behaves  like  an  ideal  fluid,  and  cannot  support  any  shear  stresses,  in 
order  to  alleviate  this  problem.  Eq.  ( 3.3 )  was  modified  to 

-/i=  (1  +wru  -  Vrt).  rtg  0.95 5  0m.  18  1) 

[,  3  /  \<PoJ 

0.045rro  /  ip  V 

=  - = —  —  (1  4-  hi)".  9  >  0.955  t)„,  (8.21 


\  3  /  J  '  ^ 

where  0m  equals  the  melting  temperature  of  the  target  material.  We  note  that  the  numerical 
simulations  by  Chen  [3]  ot  penetration  tests  involving  conical-nosed  penetrators.  and  the  cavity 
expansion  model  of  Forrestal  et  al.  [:>]  do  not  consider  thermal  softening  of  the  material. 
However.  Chen  and  Batra  [21j  have  shown  that  thermal  softening  influences  strongly  target's 
deformations. 


3. 1  Results  for  a  hemispherical-nosed  penetrator 

Figure  2  exhibits  the  computed  depth  of  penetration  both  with  and  without  the  consideration  of 
frictional  forces  at  the  target  penetrator  interface,  and  the  test  values  from  Tables  1  and  2  of 
Forrestal  et  al.  [5]  for  r0  =  2.o4mm  and  3.555  mm.  The  abscissa  equals  the  non-dimensional 
number  x  =  o(t’0:'<r0  sometimes  also  referred  to  as  the  damage  number  [27],  Higher  values  of 


- Results  without  friction 

-  Results  with  friction 

O  Experimental  values  from  Forrestal  et  al.[5] 


Fig.  2  a 


37 


d  e 

Fig.  3.  Tunnel  shapes  with  ^  =  0.12  when  the  penetrator  speed  has  been  reduced  to  a  0.7969,  b  0.6050,  and 
c  0.306  7  times  the  initial  speed  of  1 .009  km  s.  Tunnel  shapes  at  the  end  of  the  penetration  process  for  n  =  0. 1 2 
and  0.0  are  given  in  d  and  e  respectively 


A  comparison  of  these  tunnel  shapes  with  those  reported  in  Figs.  5a-d  of  Chapter  3  of 
Chens  doctoral  dissertation  [28]  suggests  that  for  the  case  ot  a  smooth  target/ penetrator 
interface  more  material  on  the  sides  of  the  interface  is  deformed  severely  as  compared  to  that  in 
the  present  case.  When  frictional. effects  are  considered,  the  size  of  the  severely  deforming  region 
ahead  of  the  penetrator  barely  changes  with  time.  This  could  be  due  to  the  retardation  of  the 
relative  movement  of  the  target  material  along  the  target,  penetrator  interface.  Thus  it  will  hinder 
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Fig.  7.  Dependence  of  penetration  depth  upon  /i  for  li)  penetrator  radius  =  3.555  mm.  striking  speed 
=  1.009 km/s  (ii)  penetrator  radius  =  2.54mm,  striking  speed  =  980m, s 


Figures  8  a,  b  depict  respectively  the  time  histories  of  the  non-dimensional  axial  resisting 
force  experienced  by  the  penetrator,  penetrator  position  and  its  speed  for  different  values  of  ft. 
The  axial  resisting  force  is  nondimensionalized  by  nr01<j0,  time  by  r 0jr0,  penetrator  speed  by  v0, 
and  the  penetrator  position  by  its  length.  Here  cr0  equals  the  yield  stress  of  the  target  material  in 
a  quasistatic  simple  compression  test.  In  the  beginning  the  resisting  force  is  minimally  affected  by 
the  value  of  ft.  However,  after  about  10%  of  the  duration  of  the  penetration  process,  the  resisting 
force  gradually  increases  to  about  10%  higher  than  that  computed  without  the  consideration  of 
the  frictional  force.  This  increase  in  the  resting  force  is  caused  by  the  tangential  traction  exerted 
by  the  deforming  target  material  on  the  penetrator  surface.  The  normal  traction  on  the 
target,  penetrator  interface  is  dominated  by  the  hydrostatic  pressure  which  is  affected  very  little 
by  the  consideration  of  frictional  forces.  As  expected,  an  increase  in  the  value  of  ft  increases  the 
axial  resisting  force  and  thus  decelerates  the  penetrator  quicker.  For  the  steady  penetration 
problem,  Chen  and  Batra  [4]  found  that  the  axial  resisting  force  depended  moderately  upon  ft. 
The  oscillations  in  the  axial  resisting  force  are  probably  due  to  high  frequency  oscillations 
resulting  from  multiple  stress  wave  reflections  from  the  top  surface  of  the  target.  In  our  work  we 
have  not  introduced  any  artificial  viscosity  to  smoothen  out  these  oscillations.  The  penetrator 
deceleration  computed  by  Chen  [3]  also  exhibited  oscillatory  behavior.  However,  the  decelera¬ 
tion  computed  by  differentiating  twice  with  respect  to  time  the  penetrator  position  was  found  to 
be  smoother  than  that  computed  directly  from  the  numerical  solution  of  the  problem.  For  the 
case  of  the  rough  target,  penetrator  interface  the  penetrator  speed  decreases  gradually  to  zero 
even  towards  the  end  of  the  penetration  process.  However,  for  smooth  interface,  the  rate  of 
decrease  of  the  penetrator  speed  drops  sharply  towards  the  end  of  the  penetration  process. 
Because  of  the  higher  deceleration  of  the  penetrator  at  larger  values  of  n,  the  total  duration  of  the 
penetration  process  decreases  as  ft  is  increased. 
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Fig.  9.  Dependence  upon  r„.  r0  of  a  [he  penetration  depth  and  b  the  total  duration  of  the  penetration  process 
for  n  =  0  and  0.12  at  initial  striking  speed  r„  of  1.009km  s.  and  for  u  =  0.0  at  r0  =  si9/n  s 

Values  of  rni  r0  close  to  0.2  imply  a  blunt  nose,  and  those  near  2  correspond  to  a  moderately  sharp 
nosed  penetrator.  In  the  plots  of  Fig.  9.  the  penetration  time  and  the  penetration  depth  are 
nondimensionalized  with  respect  to  their  values  lor  the  smooth  hemispherical  nosed  penetrator. 
For  n  =  0.  the  nose  shape  affects  the  normalized  penetration  depth  and  the  duration  of  the 


104 


R.  C.  Batra  and  X.  Chen 


periphery  essentially  equalled  zero  for  each  nose  shape  studied.  The  intense  deformations  of  the 
target  were  found  to  spread  deeper  into  the  target  for  the  blunt  nosed  penetrator  as  compared  to 
other  nose  shapes  which  agrees  with  the  presently  computed  results. 

For  many  of  the  penetration  problems  studied  herein  a  thin  layer  of  the  target  material 
adjacent  to  the  target/ penetrator  interface  melted  and  thus  underwent  severe  plastic  deforma¬ 
tions.  However,  ahead  of  the  penetrator  nose  no  thin  intensely  deforming  regions,  usually  known 
as  adiabatic  shear  bands,  were  observed.  This  is  probably  due  to  the  fact  that  no  failure  criterion 
has  been  induded  in  the  problem  formulation  and  that  the  target  has  been  taken  to  be  very  thick. 


4  Conclusions 

We  have  studied  the  dynamic  thermomechanical  deformations  of  a  very  thick  thermally 
softening  viscoplastic  target  impacted  at  normal  incidence  by  a  cylindrical  rod.  The  frictional 
force  on  the  target,  penetrator  interface  has  been  modelled  by  a  velocity  dependent  relation 
proposed  earlier  by  Chen  and  Batra  [4].  Results  for  different  values  of  the  coefficient  of  friction 
and  penetrator  nose  shapes  have  been  computed  and  presented.  The  computed  depth  of 
penetration  has  been  found  to  agree  very  well  with  the  test  values  reported  by  Forrestal  et  al.  [5]. 
For  the  same  striking  speed  the  penetration  depth  increases  rapidly  when  rnir0  is  increased  from 
0.2  to  about  0.7  and  from  1.2  to  2.0.  but  quite  slowly  for  r„,  r0  between  0.7  and  1.2.  However,  when 
frictional  effects  are  neglected,  higher  values  of  r„,r0  result  in  greater  values  of  the  penetration 
depth.  At  an  intermediate  stage  of  the  penetration  process,  more  of  the  target  material  ahead  of 
the  penetrator  and  on  its  sides  has  been  deformed  severely  for  a  blunt  nosed  penetrator  as 
compared  to  that  for  penetrators  with  other  nose  shapes.  In  each  case  a  thin  layer  of  the  target 
material  adjoining  the  target  penetrator  interface  essentially  melted. 
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